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Abstract 

The generalized Bohr Hamiltonian is applied to a description of low- lying collec- 
tive excitations in even-even isotopes of Te, Xe, Ba, Ce, Nd and Sm. The collective 
potential and inertial functions are determined by means of the Strutinsky method 
and the cranking model, respectively. A shell-dependent parametrization of the 
Nilsson potential is used. An approximate particle-number projection is performed 
in treatment of pairing correlations. The effect of coupling with the pairing vibra- 
tions is taken into account approximately when determining the inertial functions. 
The calculation does not contain any free parameter. 

1 Introduction 

For many years the Generalized Bohr Hamiltonian (GBH) 0, ^ has been still forming 
the only one model which, having possibilities to be derived microscopically (e.g. @, ^, P]), 
pretends to be able to describe collective excitations in the transitional nuclei or nuclei 
soft with respect to the quadrupole deformations. This is because the GBH is a scalar 
under rotations and, thus, possesses eigenstates of a good angular momentum whereas 
the most of other approaches deal with intrinsic states of an indefinite spin. Also, a 
rotation-vibration coupling, so important for transitional nuclei, is automatically taken 
into account in the Hamiltonian. 

The aim of the present paper is to investigate the low-lying collective states in even- 
even nuclei from the region of 50 < Z, < 82 within framework of the GBH. This 
region of the chart of nuclides is a good laboratory to study the collective model in its 
most general form as the standard approximations of it like the rotational and vibrational 
models are apparently not applicable to nuclei in question. Since the early calculation of 
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potential energy surfaces by Arseniev et al. [0 the 50 < Z, < 82 nuclei are referred 
to as gamma-soft or susceptible to deformations leading to triaxial shapes. This is why 
the Wilets-Jean model was applied already long time ago to interpret experimental data 
[Q. On the other hand, the properties of these nuclides have been explained within the 
Davydov-Filippov model treating nuclei as rigid triaxial ones (cf. 0). The triaxiality of 
the Xe and Ba isotopes has been discussed once more in the recent work of Meyer et al. 



10|. Other theoretical approaches have also been applied to the investigated region. The 



closest to our study are the applications of Kumar's dynamical deformation model to the 



Te nuclei [|Il|] and the Frankfurt general collective model to the Ba isotopes |]I2[. Among 
methods different from ours we mention the microscopic calculations |T^ of number- and 
spin-projected two-quasiparticle states with realistic interactions (MONSTER) and the 



Fermion Dynamical Symmetry Model [H, both applied to the Xe and Ba isotopes. The 
cranked Hartree-Fock-Bogolyubov calculations of the yrast bands and E2 properties have 



been performed for the Te and Xe isotopes []T5|, . Within the Interacting Boson Model 



(IBM-1) the Xe and Ba nuclei have been treated in the S0(6) limit [T^. Also, the proton- 



neutron Interacting Boson Model (IBM-2) has been applied to the Xe, Ba and Ce isotopes 



ISf . We cite here only the newer papers on broader ranges of nuclides from the discussed 
region. 

In Section ^ we recapitulate briefly the GBH recalling all necessary definitions and 
formulae. A new basis for the diagonalization of collective Hamiltonian is presented in 
Section |^. Some details of the basis construction are given in Appendix. To construct 
the collective Hamiltonian we use a standard microscopic model: nucleons in a deformed 
single particle potential interacting via the monopole pairing forces. We take the Nilsson 
potential with the parametrization of Seo [0 for a description of the single-particle 
motion. In treatment of the pairing forces we utilize an approximate projection technique 
based on the Generator Coordinate Method (GCM) for projecting the BCS wave function 
on the correct particle number |2^. Then, we calculate the collective potential by means of 
the Strutinsky macroscopic-microscopic method |^ and the inertial functions by means 
of the cranking model All the microscopic procedure is presented in Section ^ We 
know for a long time that the pairing correlations affect strongly the collective excitations. 
In early calculations for Xe and Ba isotopes we have decreased artificially the strength of 
the pairing interaction in order to simulate the couphng between the collective quadrupole 
and pairing vibrations 0. Nowadays, when ampler data are available, we know that the 
effect cannot be approximated satisfactorily in that simple way. Sakamoto p2| has tried 
to remedy description of the collective states in xenon isotopes by including in a static 
way the quadrupole pairing interactions. We do believe that it is a dynamical effect. 
Here, we approximate the effect of coupling with the pairing vibrations by calculating the 
most probable dynamical value of the pairing energy gap as a function of deformation and 
putting it to microscopic calculations of the collective inertial functions. The procedure 
is discussed in Section |[ This way we have constructed the collective Hamiltonian and 
have solved it as it stands with no free parameters to be fitted. The results for the Te, Xe, 
Ba, Ce, Nd and Sm isotopes are presented in Section Conclusions from our research 
are drawn in Section |^. 
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2 Recapitulation of the collective model 



The five dynamical variables of the model constitute a quadrupole tensor and are usually 
denoted as a^, fi = —2, —1, ... ,2. In order to separate rotations from vibrations, the 
variables a's are parametrized in terms of two deformation parameters, /3 and 7 and three 
Euler angles, 0, 6, ip {01 Q, in short), defining the orientation of the intrinsic principal 
axes with respect to the laboratory frame: 



a, 



Djo(fi)/?cos7 + ^ {D%in) + Dl_,in))(3sm^ 



(1) 



where D^^{(f),6,ilj) are the Wigner functions. Eq. (|l]) defines the intrinsic frame up to 
24 rotations forming the octahedral group O [^, |2^. In consequence, only /5 is a unique 
function of a's. The remaining intrinsic variables undergo transformations of the group O. 

The collective Hamiltonian, when expressed in terms of the intrinsic variables, is de- 
composed as 

Hcou = XihiP, 7) + totiP, 7, fi) + VconiP, 7) , 



(2) 



where Kou is the collective potential, the kinetic vibrational energy reads 



vib 



1 



2y/wr 1/3' 



-9. 



MP\l-B,,d,]-d,[P\l-B,^d^ 



+ 



P sin37 
and the rotational energy is 



sm3-fBp^dpj + -d. 



tot — 7^^ ^kl -^k 
^ k=l 



(3) 



(4) 



77' 



The mass parameters (or vibrational inertial functions) are denoted as B^^ and B^ 
They are, in general, functions of (3 and 7. Moments of inertia J7fc, (/c = 1, 2, 3), in general 
also 7)-dependent, are conventionally expressed by rotational inertial functions B^ : 



Jfc = 45fc(/3,7)/32sin2(7 -27rA;/3) for A; = 1,2, 3. 

The intrinsic components of total angular momentum are denoted as Jfc, {k 
The determinants of the vibrational and rotational tensors are: 



w 



5^7' ^ 



The volume element has to have the form 

dr = l3'^\/wr\ sin37|d/?d7di7 



B1B2B3. 



wrdTQ 



(5) 
1,2,3). 

(6) 
(7) 



in order to fulfill the hermiticity condition for TicoU- The collective electromagnetic E2, 
Ml and EO operators are determined by their (/?, 7)-dependent intrinsic components in 
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the following way: 

M(EO) =ZeR\(3,-f). 

Although the above description of the nuclear collective excitations is sometimes called 
a geometrical model it is needless at this stage to refer to the geometrical picture of nucleus. 
The model is just defined by seven functions of /3 and 7: Koih the collective potential 
and -B/3,/3, -B/3,7, -87,7, -Bi, B2, -B3, the inertial functions. To investigate electromagnetic 
properties additional functions should be defined, for instance, Qo and Q2, the quadrupole 
electric moments, Qi, Q2 and ^3, the gyromagnetic functions and i?^, the square of electric 
radius, as presented in eq. (|^) above. All of these 13 functions have to be determined 
from a microscopic theory. The geometry may conceivably be inherent in there. 



3 The basis in the collective space 



The methods used hitherto to solve the eigenvalue problem for Hamiltonian of eq. can 
be divided into two groups. The former group consists in a direct numerical solution of 
a system of partial differential equations using finite difference methods [Q, ^ The 
latter methods resolve themselves to matrix eigenvalue problems by means of expanding 
the eigenfunctions in a complete set of functions of /3, 7, 0, 6*, ip ^ |29|. Then, 

the problem of construction of an appropriate basis for a given angular momentum arises. 
The five-dimensional oscillator wavefunctions of a good angular momentum constitute a 



natural basis [El, pO, 131], p2|, p3|. 



Following the idea of Libert and Quentin we construct the basis from a complete 
set of square integrable functions of [3q = /9cos7, P2 = /5sin7 and Q in the form of 
product of gaussians and monomials in /5's, and the Wigner function in Q: 



-^^'i\(3,r{(32rDi,^{n) 



(9) 



for nQ,n2 = 0, 1, 
set of functions 



0,1, 



and M, L 



I. This basis is equivalent to the 



cosm7 
sin m7 



(10) 



with n = 0, 1, . . . and m = n, n—2, . . . , or 1. The functions ^^Lmnil^i 7' ^) form a com- 
plete set of nonorthogonal, square integrable functions with respect to measure dro. The 
basis functions (|10D depend on a parameter /i which is to be chosen further. Then we 
symmetrize appropriately the functions 99 of eq. (0) under the octahedral group O using 



a method close to that presented in ref. [34]. Details of the symmetrization procedure are 
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published elsewhere The 0-invariant functions generated from the set of eq. (|T0]) 

and forming a basis take the following form: 



Lmn 



iC=even>0 



with ranges of J, M, n and m given above. Quantum number L runs over natural even 
values within a range specific for given / and m [^. The form of functions fimKil) 
possible values of L are given in Appendix. The Euler angles dependence of \^'s is given 
through the Wigner functions: 



21 + 1 



167r2(l + 6_ 



;i2) 



The nonorthogonal basis of eq. 
instance, Kumar does 



28 



PI ) can be used to diagonahze Hamiltonian TicoU as, for 
Instead, we orthogonalize it applying Cholesky-Banachiewicz 
procedure |^6[ as is described in Appendix. In order to make the basis most effective the 
parameter fi has to be chosen appropriately. The natural choice is to fit potential V^cou 
to a paraboloid of revolution with stiffness C, to take a mean value B of all of inertial 
functions and then to use the oscillator formula: 



(f) 



1/4 



(13) 



We have tested the above choice of fi and checked that it works fairly well. The maximal 
value of the order of the polynomial /?" in eq. (|lT]) is fixed by the plateau condition. In 
our case it was n = IQ. 



4 The microscopic model 

We describe microscopically the nucleus as a system of nucleons moving in a deformed 
mean field and interacting through monopole (state- independent) pairing forces. Thus, 
the microscopic Hamiltonian reads: 

H = 4.p. + impair (14) 

with 

4.p. = Y.(^\h\u){ctc, + c\c.,) , (15) 

u>0 

^pair = - E G',E(H(^+tr3)|z/)c+cl,E(^'l(^+^^3)kV.'C.', (16) 
t=±l/2 u>0 ^ v'>Q ^ 

where t = ±1/2 is the isospin projection of neutron and proton, respectively. We use the 
Nilsson single particle Hamiltonian h with the shell dependent parametrization of the 
Is and F correction terms. Namely 

Vis = -2huoKNd-s, (17) 

Vi2 = -huJo{flNd'^ - </iAr«l^>shell), (18) 
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where 



Km = ^o[l + 8fiNi{N + 3/2)] + t^^A'/^P.^m, (19) 
m = f^oPlm (20) 

and Ps^Ni, Pi,Ni are the following integrals (evaluated numerically) of the harmonic oscil- 
lator radial wave function RNi{r): 

Pi,Ni = / Rliirydr (21) 

Jo 

fRo+a/2 

Ps,m = / R%,{rydr (22) 

jRo-a/2 

Here Rq and a stand for the nuclear radius and the surface thickness, respectively. We 



take, after ref. [T^, kq = 0.021, ki = 0.9, yUo = 0.062 and we take A = 130 as the average 
nuclear mass in the considered region. 

When solving the eigenproblem of the pairing Hamiltonian impair we replace the stan- 
dard BCS formalism by an approximate projection of the BCS wave function on a given 
particle number ||2^. This procedure is based on the gaussian overlap approximation 



(GOA) of the generator coordinate method (GCM) and gives results close to the exact 
projection. Moreover, the GCM projection is very convenient because its results are easy 
to present in terms of the usual BCS formalism. The form of BCS equations for the energy 
gap A and the Fermi energy A remains unchanged. The only effect of the particle number 
projection is that the standard BCS energy (of particles with isospin projection t) 

rBcs(A) = Y.i^ey^ - Gtvt) - ^, (23) 
is now corrected so as that the pairing energy is 

V...,(A) ^ V^ecs(A) - (e - ^><"- - -a + ^^"^'V + G^ ^ ^ ^^^^ 

Here is an eigenvalue of h and vl, u1 = \ — are the usual BCS occupation proba- 



bilities, and Ey = J (e^ — A)^ + is the quasiparticle energy. The 2Jl5Z (or A^) levels 



closest to the Fermi surface were taken when solving the BCS equation. It was shown |^ 
that when the projection on a given particle number is performed the pairing strength 
parameter should be renormalized with respect to the standard value in order to get the 
proper energy gaps. With the estimation: 

_ r go/{A-Zf/' for t = 1/2, 
^* " \ go/Z^' for t = -1/2 ^^^> 

where go = 0.285 huo, we can reasonably reproduce the experimental values of energy gaps. 
This is shown in Fig. |T] where experimental gaps deduced from measured or estimated 



38 1 nuclear masses are compared with theoretical values of neutron gap A^ 
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Fig. 1: The theoretical and experimental neutron energy gaps for isotopes of Te, Xe, Ba and Ce. Exper- 
imental points obtained from measured or estimated |]38| nuclear masses are represented by full or open 
triangles, respectively. Theoretical values are marked with lines. 



We apply the Strutinsky microscopic-macroscopic method |^ to calculate the collec- 
tive potential and we use the cranking model for determining moments of inertia, gyro- 
magnetic functions and vibrational inertial functions. The standard cranking formulae 
read [g i: 



[U„V„, - U„fV„ 



/\ |2 



'3k = ^ 2^ — ^ , ^ {T^\3k\v){vW\v) 



(26) 



(27) 



where jk and /ifc (/c = 1, 2, 3) are components of the single particle total angular momen- 
tum and magnetic moment, respectively, and 



Jo {u^v„' + Ui,>Vu)^ dh dh 



u,u'>0 



1 

+ - V — 

^ iy>0 



' dX , ,dh , , dX , ,dh , , 



dq' 



dX dX 



dq dq' \ dq dq' 



dq' dq 



(28) 



for q, q' = P or 7, where /(/5) = 1, /(7) = p. To be complete we quote at the end 
formulae for the electric moments: 



Qo= ZeJ^Y.^vl{u\\{l-nyY,o\u) 
y ^ u>o ^ 



(29) 



7 



Q2 



/ 1 f^TT 1 t'^ 



u>0 



V2' 



(30) 
(31) 



5 Effect of the collective pairing vibrations 



It is suspected for a long time ^ that the pairing vibrations are coupled to 

the collective quadrupole motion. Here we do not construct a complete "quadrupole + 
pairing" collective model which would have nine degrees of freedom and could not be easy 
to solve. We just take approximately into account the effect of the pairing vibrations on 
the quadrupole excitations in question. 

The collective pairing Hamiltonian has the following structure [O, M : 



n 



d A/wA(A)r$(A) d 



pair 



+ 



2y'u;A(A)r$(A)^A 5aa(A) 9A 2^^) 

+Ap,i,(A)A^+ypai,(A) 



(32) 



where J\ = —id/d^ is the particle number excess operator and $ is the gauge angle. 
Expressions for the functions i?AA(A), J7$(A), W/\{A), r$(A) andylpair(A) obtained from 
GCM and GOA are given in ref. pf. 




0.00 0.12 0.24 0.36 0.48 0.60 0.72 



Fig. 2: The neutron zero-point pairing vibration at deformation /? = 0.15, 7 = 60° for neutron number 
N = 74. The equihbrium value of the energy gap Abcs ~ 0.20?iti;o while the most probable value 
Ao w 0.12?iwo- 
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0.63 




Fig. 3: Ratio of the neutron energy gaps Aq/Abcs as a function of /3 and 7 for the neutron number 
N — 74. The minimal and maximal values of the ratio are marked on the surface. 



The ground-state energy, Eq, of the pairing vibrations is to be found from the following 
eigenvalue problem: 

7ipair^o(A)=^o*o(A), A^^o(A) = 0, (33) 

which we solve at each deformation point 7) for a given number of neutrons or 
a given number Z of protons, respectively. This way we find i?o(/3, 7) and \E'o(A;/5,7) 
separately for neutrons and protons in a given nucleus. Then, we replace Abcs of the 
BCS (such that 1/pair(A) = min) in all collective functions (eqs (H), (|^), (||) and (|2|)) 
by the most probable value of the energy gap Aq, i.e. for which the probability 

P(A) = l^or^Ar^ (34) 

of finding a given value of A in the ground state is maximal (Fig. Since -Baa decreases 
rapidly with the increase of A, the values of Aq are systematically lower than the equi- 
librium energy gaps for all values of /? and 7 (cf. [0). An example of such a behaviour 
is demonstrated in Fig. |^. 



6 Results 

The collective model in version presented in Section |^ with the Hamiltonian and elec- 
tromagnetic multipole operators constructed in the framework of microscopic model dis- 
cussed in Section ^ has been applied for a description of the collective states in even-even 
isotopes of Te, Xe, Ba, Ce, Nd and Sm. The effect of coupling with the pairing vibrations 
as reported in Section ^ has been taken approximately into account. The calculations have 
been performed with no free parameters to be fitted. The results are discussed below. 
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6.1 The energy spectra 

The calculated energy levels in the six isotope chains in question are shown in Figs ^ ^, 
I, 0, H and H, respectively and compared with the experimental spectra which are taken 
from MM unless another reference is quoted. 
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Fig. 4: Calculated (th) and experimental (exp) energy levels in the even-even samarium isotopes 
i32-i42gj-|^^ The theoretical levels are marked with spin (parity of all of states is positive), while the 
experimental ones — with spin and parity. Experimental data are taken from . 
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Fig. 5: The same as in Fig. but for neodymium isotopes 
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Fig. 6: The same as in Fig. Wbut for cerium isotopes Ce 
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Fig. 7: The same as in Fig. ^ but for barium isotopes ^^^^^^^Ba. 

The energy spectra in all of the isotopes of the two heaviest elements in question, Nd 
and Sm, are reproduced fairly well. The yrast states 2"*", 4"*", 6+ and 8"*" and the second 
2+ state are in more or less correct positions, although the yrast bands are little bit too 
stretched (Figs ^ and |^). Similar results we have got for three lightest cerium isotopes. 
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Fig. 8: The same as in Fig. but for xenon isotopes 
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Fig. 9: The same as in Fig. W but for teUurium isotopes 
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Te. 



^^"^Ce, ^^^Ce and ^^^Ce while the first 2"^ state becomes too high for heavier isotopes. In 
the case of two heaviest ones, ^^^Ce and ^^^Ce, also two-phonon triplet of states 0+, 2+, 4+ 
is visibly too high (Fig. H). The spectra only of two or three lightest isotopes of barium, 
xenon and tellurium, ^^^Ba, ^^^Ba, ^^^Ba, ^^'^Xe, ^^°Xe, ^^^Te, ^^^Te, are of a good energy 
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Fig. 10: The same as in Fig. ^ but for light barium isotopes ^^"^ ^^°Ba. Experimental data for -'^-'^^Ba are 
taken from iBsll. 



scale. For the heavier isotopes the calculation gives too stretched energy spectra (Figs 
0, H and P). The theoretical energy levels for very light barium isotopes ^^"^"^^^Ba are 

The results for ^^^Ba are compared with very recent 
Here the energy scale of the theoretical results for 



plotted additionally in Fig ITO 



experimental data from Ref. ^ 



^^^Ba and ^^°Ba seems also to be too stretched with respect to the experimental tendency. 
Probably the same occurs for unknown yet levels in ^^^Ba and ^^^Ba which have recently 



been investigated for the first time 



6.2 The electromagnetic properties 

The results of calculation of the reduced transition probability from the first excited 2^ 
state to the ground state, -B(E2; 2+ 0+) are compared with experimental data in Fig. |11 
It is seen that the 5(E2) values for isotopes of Ce, Nd and Sm are reproduced quite well. 
The calculated probabilities for most cases of the Xe and Ba isotopes are only a little 
too large and a tendency to decrease probabilities with increase the neutron number is 
reproduced correctly. An agreement of the theory with experiment is, in general, better 
for heavier Xe and Ba isotopes and worse for lighter ones. A similar observation can be 
made for tellurium isotopes. However, the theory gives i?(E2) values plainly too big in 
this case although the tendency of changes from isotope to isotope is again correct. 

Theoretical values of the spectroscopic quadrupole moment of the first 2"'" state are 



shown in Fig. 0. The experimental data are scanty for this quantity and it is difficult 
to estimate of agreement between calculations and measurements. In cases of ^^°Ba and 
^^^Ba it is very good. An opposite sign of the quadrupole moment is predicted for ^^^Ba. 
For the Te isotopes the agreement is rather poor. 
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Fig. 11: Calculated and experimental |Q reduced transition probabilities B(E2; 2+ 0+) (in c^b^) from 
the first excited 2+ state to the ground state. 
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Fig. 12: Calculated and experimental [Q spectroscopic electric quadrupole moments (in eb) of the first 
excited 2+ state. 



The calculated magnetic dipole moments agree with the experimental values in all 
of known cases of the Xe, Ba and Nd isotopes. The agreement between the theory and 
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experiment for the tellurium isotopes is again poorer. This is shown in Fig. 



|T3| 
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A 



Fig. 13: Calculated and experimental magnetic dipole moments (in /in) of the first excited 2+ state. 



Fig. |Tj shows the results of calculation of isotopic shifts for the charge root mean 
square radii. The experimental values in all of cases available for the Xe and Ba isotopes 



[J7| are reproduced very well. 
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The knowledge of experimental data for other electromagnetic properties like transition 
probabilities from higher excited states, branching and mixing ratios is still rather poor 
in the investigated region of nuclear chart and, therefore, we do not discuss them here. 



7 Conclusions 

The low-lying quadrupole collective states of 48 even-even nuclei from the region 50 < 
Z, N < 82 of the chart of nuclides are investigated within the framework of the gener- 
alized Bohr Hamiltonian. The collective potential and inertial functions are determined 
microscopically. There are no free parameters in the calculation. Taking the coupling 
of the quadrupole and pairing vibrations into account brings the energy levels down to 
the scale comparable with that of experimental one. In general, the GBH works better 
in the case of elements with larger Z, namely, Ce, Nd and Sm than for Xe, Ba and, 
especially, Te which has only two protons outside the closed shell Z = 50. Energies are 
better reproduced by the calculation for isotopes with lower number of neutrons. For 
those with neutron number = 78, 80 the energy levels rule, too high. On 

the contrary, electromagnetic properties seem to be better reproduced just in the case of 
heavier isotopes. The calculated reduced probability of the 2+ — O"*" transition, if not 
reproduced well, has a tendency to be too large, especially in cases of lighter isotopes for 
which energy levels lie in more or less correct position. 

It does not make sense to discuss how good or how bad are our results obtained with 
no free parameters in comparison to those of other approaches in which some parameters 
are usually fitted to the data. 

Finally, we may say that though the collective model based on the generalized Bohr 
Hamiltonian gives qualitatively reasonable results for nuclei from the major shell 50 < 
Z, N < 82, it is not able to reproduce all of collective characteristics simultaneously in 
detail. We think that dynamical effects of coupling with pairing vibrations can improve 
the description of collective states in question. 

Appendix. 

Symmetrization under the O group 

All of transformations which change names and senses of the intrinsic axes and form 
the octahedral group O can be expressed through the three well known Bohr rotations, 
Ri, i?2, Rs (see refs [|| |^, H, ^ for their definition). Therefore, it is sufficient to sym- 
metrize functions 7, fi) of eq. ( ]TD| ) with respect to Ri, {i = 1, 2, 3). Symmetriza- 
tion under Ri and R2 gives 



il + R,)il + R, + Rl + Rl)^ZniP.7,^) oc = e-'^^V2/3X^^(7)$k(^^) (35) 



where 




cosm7 for L/2 even, 
sinm7 for L/2 odd 



(36) 
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and ^Mi^ is given by eq. (|T2]). The additional symmetrization under R3 leads to the 
functions: 

(1 + i?3 + RD^l'LiP, 7, ^) = ^ilniP, 7, ^) = e-^^'/'(3- Y: fLKil)KlKi^) (37) 

_fsr =even>0 

where 

fimAl) = CLmKPKmil)- (38) 

Coefficients cj;^rnK equal to 



cimK = + 2(-l)^/^[(l + 5,,)(l + <5^o)]-'/'4l(^/2) 



2cos(2m7r/3) for (iT + L)/2 even, 
flLmK = 2sin(2m7r/3) for K/2 even, L/2 odd, (39) 
-2sin(2m7r/3) for fs:/2 odd, L/2 even. 

Function d^Kii^) is related to the Wigner function in a usual way: Dj^i^{(j),9,ip) = 

Since i?3 does not commute with Ri and R2 the set of functions of eq. (|37[) is linear 
dependent. Selection of linear independent functions from the set ( P7| ) leads to restrictions 
on values of quantum the number L for given / and m. The number L takes even integer 
values ranging from 25/ odd = 2(1 — 5/ even) to L^^.^^{m), where 5/ even is equal to 1 for an 
even / and for an odd one, and 

2 ■ [7/12] + 2Aj - 25/even for m = 0, 
2 ■ [1/6] for m > 0, /mod3 = 1, 



max 



with 



2 ■ [7/6] — 25/ even for m > 0, mmod3 = 0, J mod 3 = 2, 

or mmod3 7^ 0, / mod 3 = 0, 

2 ■ [//6] + 25/odd for m > 0, mmod3 = 0, Jmod3 = 0, 

or m mod 3 7^ 0, J mod 3 = 2, 

(40) 

_fl for /mod 12 = 0, 4, 6, 8, 9, 10, , . 

^"10 for /mod 12 = 1, 2, 3, 5, 7, 11. ^ ^ 



A rational number N/M in square brackets, [N/M], stands for its integer part. 

We have to consider also the boundary conditions in 7 (it is for 7 = kn/S). The 
combinations of the linear independent functions ^E^i^n which fulfill these conditions and, 
thus, form a (nonorthogonal) basis in the domain of the Bohr Hamiltonian, are (cf. 



^iZniP,l,^)=e->'^'/'p- E fLMKiK (42) 

K=even>0 

where fl^K the following combinations of functions fl^K of ^l- dH) for all possible 
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values of K: 



m < 6 



fa 



OOK 



fooK 

fo2K — fo2K + 2^00K 
fL3K = flSK + 25lo/o1A' 

fU = fU - {-^f''fL2K for 
ILk = ILk - {-^f'^fU for 



mmodS 7^ 

fLmK = fLmK ~ fim-QK fo^ 

mmodS = 

flmK = /imi^: ~ ^ L mod 4. ofim-QK fo^ 



L 
L 



(0),2,6,10 
(0),2,4,.., 
(0),2,4,.., 



,4-[(LLx(3) + 2)/4]-2 
(2) 



ji 

' max 
' -^max(l 



m > 6 
L = (0),2,4, 
L = (0),2,4, 



^max("^) 



,^max("^) 



where the lowest value of L in parentheses is omitted for odd J's. 

The orthonormal (with respect to the measure dro) set of functions reads: 



where 



and 



L'm' 



(43) 

(44) 
(45) 

(46) 



Lower diagonal matrices S^"^\jj?j and -E'-^-* are obtained through Cholesky-Banachiewicz 
method. The running index m' in eq. (^) has the same parity p as that of m. So, we 
have two matrices E^^p\ one for each parity, for a given I. On the other hand, matrices 
S^"^\fi) depend on m itself. They depend also on parameter ft. However, it turns out 
that 

RnM /i) = fJ^'^'R^M; /i=l) (47) 

Hence, it is sufficient to perform the orthonormalization only once. 

Hamiltonian T^cou is hermitian with the measure dr of eq. (|^) . The basis orthonormal 
with respect to this measure is finally of the following form: 



:pIM 
Lmn 



(4J 
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